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Abstract 

We prove the finite time blow-up for C 1 solutions to the Euler- 
Poisson equations in R n , n > 1, with/without background density for 
initial data satisfying suitable conditions. We also find a sufficient 
condition for the initial data such that C 3 solution breaks down in 
finite time for the compressible Euler equations for polytropic gas 
flows. 
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1 The Euler-Poisson equations 

We are concerned on the following (pressureless) Euler-Poisson equations in 
W n , n>l. 

' d t p + div (pv ) = 0, 
d t (pv) + div (pu <g> u) = —kpV(f), 
A0 = p, ( ' } 

k v(x, 0) = v (x), p(x, 0) = p (x). 

The unknowns are the velocity field v — (v 1 , ■ ■ ■ , v n ) = v(x, t), mass density 
p = p(x,t) > 0, and the potential function of the force field = <p(x,t). We 
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consider the case k > 0, which represents that the force is attractive. Phys- 
ically the system (11.11) represents the large scale dynamics stars or cluster 
of stars, regarded as particles consisting of a gas, interacting through grav- 
itation. The global regularity/finite time blow-up problem from the view 
point of the 'critical threshold phenomena' have been studied extensively by 
Tadmor and his collaborators for the one dimensional Euler-Poisson equa- 
tion (with/without pressure term) (see e.g. OH], [12] and references therein). 
In the multi-dimensional problem, the finite time blow-up for the attractive 
force case with pressure in the spherical symmetry for compactly supported 
p(x, t) in IR 3 was obtained in [7] , and for the repulsive case with similar ge- 
ometry the blow-up was deduced in [ID] . On the other hand, the global 
regularity for some large class of initial data near a steady state is obtained 
in [5]. In [9] the stability type of result was obtained with inclusion of the 
pressure. As for the study for the equations (11.11) in the general setting, 
which is our main concern in this paper, the main difficulty was the forcing 
term fcV0, which is of nonlocal nature, and resembles the notorious pressure 
term in the 3D incompressible Euler equations. This feature was emphasized 
in [8], and was the main motivation for studying a model problem, called 're- 
stricted Euler-Poisson equations', where the nonlocal forcing term is replaced 
by a local one. The finite time blow-up for the restricted Euler equations was 
shown via spectral dynamics in [8]. Our aim in this paper is to show that the 
finite time blow-up is unavoidable for the original Euler-Poisson equations 
under suitable conditions for the initial data of large class. At least for our 
blow-up problem the nonlocality can be bypassed simply by decomposing 
the evolution equations of the velocity gradient matrix into the symmetric 
part (deformation tensor equation) and the antisymmetric art(vorticity equa- 
tion), and observing that the vanishing property of the vorticity is preserved 
along the particle trajectories, concentrating on the evolution of the den- 
sity, not the velocity gradients. Let us define the vorticity matrix Q = (Oy) 
for the n dimensional vector field v = (f 1 , ■ ■■ ,v n ) as := — djV 1 , 
i,j — 1, • • • , n. In the 1-D case we set Q = 0. We state our main theorem 
for the system (11. ip . 

Theorem 1.1 Let n > 1 in the Euler-Poisson equations U.l\) . We suppose 
the initial data (v , p ) satisfies the following condition, 



St = {a e W 1 1 ft (a) = 0, -divt» (a) > y ^ v 1 > 0} ^ 0. 

Then, the C 1 — regularity of local classical solution with initial data t>o,Po 
cannot persist arbitrarily long time. 

Remark 1.1 We observe that in the 1-D case and in the n-D spherically 
symmetric case the vorticity vanishing condition, Qq(cl) = is redundant. 
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Remark 1.2 As will be seen on the proof in the next section the proof 
does not depend on the specific domain, and the same results hold true for 
the bounded domain with smooth boundary or the periodic domain. 



Next, we consider the following Euler-Poisson equations with constant back- 
ground density, p. 

' d t p + div (pv) = 0, 
d t (pv) + div (pu ® u) = — kpV(j), 
A0 = p - p, 
K v(x, 0) = v (x), p(x, 0) = p (x). 

The inclusion of constant density in the Poisson part of the equation has 
more physically realistic meaning in many cases. For example, if we consider 
the above system in the periodic unit box T n = W 1 fL n , then the natural 
choice of p is 

p= p(x,t)dx, 



which is constant in time thanks to the first equation of (11.21) . 

Theorem 1.2 Let n > 1 in the Euler-Poisson equations ( fi.ffj) . We suppose 
the initial data (fo,Po) satisfies the following condition, 

S2 '■= {the set of all points a G R n satisfying U.3\) - fL~4\ ) below} ^ 

where 

Q (a) = 0, divv (a) < 0, p {a) > + ^ ) p, (1.3) 

and 

'k( Po (a)-D [A(p (a) - I f - 3p 2 ] 3fcp 4 



(divvo(a-)) > max 



6p (a) 2 ' 4pg(a)(p (a) - f 

(1.4) 

Then, the C 1 — regularity of local classical solution with initial data v 0, po 
cannot persist arbitrarily long time. 

2 Proof of Theorems 1.1 and 1.2 

Proof of Theorem 1.1 We write the first two equations of (II .ip in a more 
convenient form, 

d t p+ (v ■ V)p = -pdivv, (2.1) 

and 

d t v + (v ■ V)u = -W(f) (2.2) 
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respectively. Taking a partial derivative of (12.21) . we obtain the matrix equa- 
tion, 

d t V + (v -V)y + V 2 = -fc$, (2.3) 

where we set V = (<9jt> J ), and $ = (didj<f)). The symmetric and the antisym- 
metric parts of (12. 3p are 

—V = -V 2 -A 2 - fc$, (2.4) 

and 

—A = -VA - AD (2.5) 
Dt K ' 

respectively, where we set V = \{V + V T ), A = |(V — V T )(= and 
denoted = 9 t + (u • V). We consider evolution along the particle tra- 
jectory {X(a,t)}, defined by a solution of the following ordinary differential 
equations, 

= v(X(a,t),t), X(a,0) = a e S 1 , 

where v(x,t) is a classical solution of the system (11.11) . From the equation 
(12.51) we have immediately 

— \A\ < 2\V\\A\, 



where we used the matrix norm, \M\ := y Y^ij=i ^ij- By Gronwall's lemma 
we obtain 



\A(X{a,t),t)\ < |A (a)|exp 



2 f \V(X(a,r),r)\dr 
Jo 



(2.6) 



Since A (a) = |fi (a) = for a G <Si, we have A(X(a, £),£) = along the 
particle trajectory {X(a, t)}. Hence, taking trace of (12. 4p . we find 

§- t (Tr(S)) = -\V\ 2 -kA<f> 

along the trajectories {X(a, t)}, which taking into account of the fact divt> = 
Tr(V) and the Poisson part of (11.11) . leads to 



R-(dwv) = -\V\ 2 -kp. (2.7) 



Dt 

From (12.11) . using (12.71) . we have 

D 2 p D , (Dp\ , D 



m - Dt (p divv )--{-m) divv - p Dt {divv) 

= p(divv) 2 + p\V\ 2 + kp 2 

> kp 2 . (2.8) 
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Integrating this inequality along the particle trajectory once, we have 

^(X(a,t),t)>f^\(t = 0) + kj\ 2 (X(a,s),s)ds, (2.9) 
which implies that 



if 

Dp 



^(X(a,t),t)>0, (2.10) 



1 (t = 0) = -p (a) dwv (a) > 0. (2.11) 

Note that this is the case for a & Si. Hence, multiplying ( 12. 8ft by without 
changing direction of the inequality, we find 

2Dt \DtJ ~ 3 Dt ' V ; 

and integrating this along the particle trajectory again, we have 

Dp\ 2 ^ 2k 3 /Dp _ \ 2 2A; 3 



if 



> |p 3 , (2.13) 



±£(t = 0)J - y pg(a) = p§(a)(divi*,(a)) 2 - y p^(a) > 0, (2.14) 
which satisfied in our case a G S%. The reduced inequality 

5e > Ai 

~ V 3 P 

from ( I2.13P can be immediately solved to yield 

p(X(a, t),t)> Po{a) , >9 , (2.15) 

fcpo(a) y. 
12 1 



which shows that p(X(a, t),t) blows-up at a time no latter than 
if the C 1 solution (v(x, t), p(x, t)) persists until that time. □ 
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kpo(a) ' 
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Proof of Theorem 1.2 In this case the proof is the same as that of 
Theorem 1.1 up to the parts before (12.81) . Here, instead of (12.81) we have 

° 2 P D (A' \ (°P\ rV D <rV \ 

- {p div v) — — I div v — p — (div v) 



Setting 



we find that 



Dt 2 DV r ' \Dt J Dt 

= p(divv) 2 + p\V\ 2 + kp(p - p) 

> kp{p-p) = k[p- P ^) (2.16) 
D 2 ,„ 2 hp 2 

Integrating this inequality along the particle trajectory once, we have 

^(X(a, t),t)> {^\ (t = 0)-^ + kJ* 6 2 (X(a, s), s)ds 

ko 2 t r l 

= -p (a)dwv (a) - + k / 6 2 (X(a, s), s)ds. (2.18) 

4 Jo 

This implies that 

^-(X(a,t),t)>0, (2.19) 

if 

4p (a) divfo(a) 
kp 2 

Note that for a £ 52 the first condition is satisfied, and we will concentrated 
on the time interval [0, t] defined by the second inequality in (12.201) . Hence, 
multiplying ( 12.171) by ^| without changing direction of the inequality, we 
find 

1D_ (DOV kDP _ kf_D9_ 

2Dt\Dt ) ~ 3 Dt 4 Dt ' 1 ' ' 

and integrating along the particle trajectory again, we have 

D0\ 2 2k „, kp 2 n (D6 , V 2k , fep 2 , . , 

> -^-o 3 --^-e+(—(t = o)) - —el{a) + -Z-9 (a) 



div„ (a) < 0, and < t < -5?^. (2.20) 



Dt J ~ 3 2 7 3 uw 2 

= ^ -^9 + p 2 (a)(divv (a)) 2 - 2 ±6l{a) + ^9 (a) 

k 

3' 



k 

> -6 3 , (2.22) 
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if 

pl(a)(divv (a)) 2 - Y 6l(a) + -^6 (a) > 0, (2.23) 

and 

> (2.24) 
Indeed, from (12.191) . we see that (I2.24p implies that 

WWM).«)>*g, (jr(M)|t) 
3 2 

as long as the classical solution is well-defined. The conditions (12.231) and 
(I2.24p are satisfied, in turn, if 

( , < ^ ~|) [4(po(a) ~ f ) 2 ~ 3p 2 ] 

(diWa)) > ^ , (2.25) 



and 



Po(a)> [l + Jhp (2.26) 



2 V 2 

respectively, which follows from the fact a G £2- The inequality 

- V 3 

can be solved to yield 

9(X(a,t),t) > 9 ° {a) x2 , 12.27) 

12 ^ 

which shows that p(X(a, t),t) = 9(X(a, t), t) + | blows-up at a time no latter 



than = y ke {a) > ^ ^he ^ som tion (v(x,t), p(x,t)) persists until that time. 
Comparing this with the second inequality of (12.201) . in order to get the finite 
time singularity, it suffices to have 



4p (a) div^o(a) / 12 



fcp 2 y ^6*0 (a 

which is true, if 

(div^(a)) 2 > g , (2.28) 

4 Po(«)(Po(a) - f) 

Both of the conditions (I2.25P and (I2.28P are satisfied for a G 1S2. □ 
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3 Remarks on the Polytropic Euler equations 



In this section we show that a similar argument as in the above section can 
be used to show a finite time blow-up result for the following Euler equations 
for polytropic gas flows in a domain in R n . 

' d t p + div (pv) = 0, 

dt(pv) + div (pv <g> v) = — Vp, 
< d t S + (v-V)S = 0, (3.1) 
p = p(p,S)=p^e s , ( 7 >1) 
k v(x, 0) = v (x), p(x, 0) = po(x), S(x, 0) = S Q . 

In the above p = p(p, S) is the pressure, and S is the entrophy. For both 
of the physical and mathematical backgrounds of the system (13.11) including 
the local well-posedness of classical solutions we refer the monograph [6] and 
the review articles [H [2] . For the blow-up problem of the system (13.11) for 
n = 3, in particular, there is already well-known result due to Sideris in 
[TT] (see also [H [2] and the references therein for the other related results on 
this problem). In this section we just would like to show another aspect of 
the blow-up problem of the Euler system, represented in Theorem 3.1 below. 
As in the previous section we rewrite the system (13 .ip in a more convenient 
form, 

r d t p + (v ■ V)p = — p div v, 

d t v + (v ■ V)v = -Vp, p = e s p 1 
* d t S+(v-V)S = 0, 

k v(x, 0) = v (x), p(x, 0) = p {x), S(x, 0) = S . 

We use the matrix notations, V{j = \(diV^ + djV 1 ), = \(diV^ — djv 1 ), and 
Pij = didjp as previously. We also denote {(Xk(x, t), e&(a;, t))} for the pairs 
of eigenvalue and the normalized eigenvector of the deformation tensor T>. 
In this section we shall prove the following theorem. 

Theorem 3.1 Let 7 > 2 be given. We suppose the initial data (vo,po,So) 
for AS.fy) satisfies the following condition, 

5 3 := {the set of all points a EW 1 satisfying Ii3.3\) - I{375\) below} 7^ 

where 

fio(o) = 0, (3.3) 
D k S (a) = OVA; G {1,2} andD m p (a) = OVm G {0,1,2}, (3.4) 

and 

3j G {1, • • • , n} such that Xj(a, 0) < (3.5) 

Then, the C 3 — regularity of solution to the system K3.2\) cannot persist arbi- 
trary long time. 
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Proof Taking the first and the second partial derivatives of the first equation 
of (Q, we find that 



D 
Dt 



djp = —(djV ■ V)p — djpdivv — pdivdjV, 



(3.6) 



and 
D 



Dt d jd k p = -(djd k v ■ V)p - (djv ■ V)d k p - (d k v ■ V)djp - djd k p div v 



—djp div d k v — p div djd k v — d k p div djV 



(3.7) 



respectively. Combining (13.61) and (13.71) with the first equation of (13.21) . we 
find that 

^(p + \Dp\ + \D 2 p\) < 2(\Dv\ + \D 2 v\ + \D 3 v\)(p + \Dp\ + \D 2 p\), 

which, after integration along the particle trajectories, provides us with 
(p+\Dp\ + \D 2 p\)(X(a,t),t) < 
< {po + \D Po \ + |£) 2 po|)(a)exp 



2 / (\Dv\ + \D 2 v\ + \D 6 v\)(X(a,s),s)ds 
o 

(3. 



Similarly, from the third equation of (13.21) we obtain the estimate, 

(\DS\ + \D 2 S\)(X(a,t),t) < 

<{\DS \ + \D 2 S \){a)exp 2 f (\Dv\ + \D 2 v\)(X(a, s), s)ds 

.. Jo 



(3.9) 



The symmetric and the antisymmetric parts of the velocity gradient matrix 
equations, obtained from the second equation of (13.21) are 



D 
Dt 



-V = -V 2 -A 2 -P, 



and 



—A = -VA - AV 
Dt 



respectively. In particular, from (13.111) we obtain easily that 



\A(X(a,t),t)\ < |A (a)|exp 



2 / \Dv(X(a,s),s)\ds 
o 



(3.10) 
(3.11) 

(3.12) 



We observe the pointwise estimate for the Hessian of the pressure, 

\P\ < Ce s (\DS\ + \D 2 S\)(p + \Dp\ + \D 2 p\), (3.13) 
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where C = C(p), and limsup p _^ |C(p)| < oo thanks to the assumption 7 > 2. 
The estimates (1319) . (EE5j) . (13TT2]) and fl37T3|) shows that 

A(X(a, t), t) = P(X(a, *),*) = (3.14) 

if a G 1S3, which we suppose from now on. Hence, along the particle trajec- 
tories {X(a,t)}, the equation ( 13.101) reduces to 

—V=-V 2 . (3.15) 
Dt y ' 

Let (Xj,ej) be an eigenvalue- normalized eigenvector pair of V such that 
Xj(a,0) < 0. Then, taking operation on T>ej = Xjej, we have 

D _\ _£>e, DXa , De 7 - 

.s? 2> ; ej + ^ = cf e ' + ^ 

Using (I3.15p . this reduces to 

-^"st'-Dh^D?- (3 ' 16) 

Multiplying on the both sides of (I3.16p . we find 

DA,- 



-A, (3.17) 



which can be solved along the particle trajectories {X(a,t)} to yield 
and shows that 

1 

A 3 (X(a, t), £) — > — 00 as £ — > £* 



Aj(a, 0) 
□ 

Acknowledgements 

This research was done while the author was visiting University of Nice. 
The author would like to thank to Y. Brenier for his hospitality and for 
stimulating discussions with useful suggestions. He would like to thank also 
to E. Tadmor for informing him of the recent results on the subject with 
references. This research was supported partially by KRF Grant (MOEHRD, 
Basic Research Promotion Fund). 



10 



References 



[I] G-Q Chen and D. Wang, The Cauchy problem for the Euler equations 
for Compressible Fluids, Handbook of Mathematical Fluid Dynamics, 1, 
(2002), Elsevier Science Ltd. 

[2] D. Christodoulou, The Euler equations of compressible fluid flow, Bull. 
A.M.S., 44, no. 4, (2007), pp. 581-602. 

[3] S. Engelberg, Formation of singularities in the Euler and Euler-P oisson 
equations, Phys. D, 98, no. 1, (1996), pp. 67-74. 

[4] S. Engelberg, H. Liu and E. Tadmor, Critical thresholds in Euler-Poisson 
equations, Indiana Univ. Math. J., 50, (2001), Special Issue, pp. 109-157. 

[5] Y. Guo, Smooth irrotational flows in the large to the Euler-Poisson system 
in R 3+1 ; Comm. Math. Phys. 195, no. 2, (1998), pp. 249-265. 

[6] A. Majda, Compressible Fluid Flows and Systems of Conservation Laws 
in Several Space Variables, Springer- Verlag, New York 1984. 

[7] T. Makino and B. Perthame, Sur les solutions a symetrie spherique de 
Venation d'Euler-Poisson pour Involution d'etoiles gazeuses, Japan J. 
Appl. Math., 7, (1990), pp. 165-170. 

[8] H. Liu and E. Tadmor, Spectral dynamics of the velocity gradient field in 
restricted flows, Comm. Math. Phys. 228, no. 3, (2002), pp. 435-466. 

[9] T. Luo and J. Smoller , Rotating Fluids with Self- Gravitation in Bounded 
Domains, Arch. Rational. Mech. Anal., 173, (2004), pp. 345-377. 

[10] B. Perthame, Nonexistence of global solutions to Euler-Poisson equa- 
tions for repulsive forces, Japan J. Appl. Math. 7, no. 2, (1990), pp. 
363-367. 

[II] T. C. Sideris, Formation of singularities in compressible fluids, Comm. 
Math. Phys., 101, (1985), pp. 475-485. 

[12] E. Tadmor and D. Wei, On the global regularity of sub-critical Euler- 
Poisson equations with pressure, J. European Math. Society, to appear. 



11 



